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Abstract. We study the bulk viscosity in neutron star cores due to modified Urea processes involving nucleons, 
electrons and muons and analyze its reduction by singlet-state or triplet-state superfluidity of nucleons. In com- 
bination with the results of our previous paper on the bulk viscosity due to direct Urea processes, a realistic 
description of the bulk viscosity in superfluid neutron star cores is obtained. Switching off direct Urea processes 
with decreasing density in a nonsuperfluid matter lowers the bulk viscosity by 3-5 orders of magnitude. The 
presence of muons opens additional source of bulk viscosity due to muon Urea processes and lowers the threshold 
density of the electron direct Urea process. The superfluidity may strongly reduce the bulk viscosity and affect 
thus damping of neutron star vibrations. 
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1. Introduction 

The bulk viscosity of neutron star matter determines 
damping of neutron star pulsations which could be excited 
during neutron star formation or triggered by instabilities 
developed during neutron star evolution (e.g., Cutler et al. 
|199(]| ). The bulk viscosity can damp gravitational radiation 
driven instabilities in rotating neutron stars and therefore 
influence the maximum rotation frequency (e.g., Lindblom 
1994 Zdunik 1996, Andersson & Kokkotas |00^. 



where I is either electron or muon and is an associated 
neutrino. Direct Urea processes are allowed (Lattimer et 
al. 1991) for some realistic equations of state at densi- 



In this paper, we continue to study the bulk viscosity 
in neutron star cores (the layers of density p ^ 1.5 x 10^^ 
g cm~'^) assuming that stellar matter consists of neutrons 
n, protons electrons e, and muons /i and taking into ac- 
count superfluidity of nucleons. Protons and neutrons con- 
stitute strongly interacting Fermi liquids while electrons 
and muons constitute almost ideal, strongly degenerate 
gases. Electrons are ultrarelativistic; muons may be non- 
relativistic or moderately relativistic, depending on den- 
sity. The bulk viscosity in this matter is mainly produced 
by neutrino processes of Urea type. These processes are 
known to be divided into most powerful direct Urea pro- 
cesses and much weaker modified Urea processes. 

A direct Urea process consists of two reactions, direct 
and inverse one. 



ties higher than certain threshold densities (of several po, 
where po — 2.8 x 10^*^ g cm^"^ is the standard nuclear mat- 
ter density). Thus, they may be open in the inner cores of 
rather massive neutron stars producing large bulk viscos- 
ity. The threshold density for the muon process is always 
higher than for the electron one. 

However, for many equations of state, direct Urea pro- 
cesses are forbidden by momentum conservation at any 
density in a neutron star core. Moreover, at p ^ 3 po they 
are prohibited for the majority of equations of state. Thus, 
they do not operate in the low and medium-mass neutron 
stars and in the outer cores of almost all neutron star mod- 
els. If so, the bulk viscosity is determined by the reactions 
of modified Urea processes 

n + N-^p + N + l + Di, p + N + l~^n + N + vi, (2) 

where N is an additional nucleon required to ensure mo- 
mentum conservation. For instance, in the npe matter one 
has two modified Urea processes, with N — n or p, which 
will be referred to as the neutron and proton branches of 
the modified Urea pr ocess (e.g., Friman & Maxwell 1979, 



p + I + i^i, p + l 



(1) 
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Yakovlev & Levenfish |l995[ ). In the npep matter, we have 
four processes corresponding to N — n ot p, and I = eor p; 
they will be labeled as Nl. The rates of the modified Urea 
processes are 3-5 orders of magnitude lower than those 
of the direct Urea processes. The modified Urea processes 
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either have no density threshold (as the neutron branch in 
the npe matter) or have much lower thresholds than the 
direct Urea. 

Sawyer (|1989D and Haensel & SchaefFer (|1992|) ana- 



lyzed the bulk viscosity of non-superfluid npe matter pro- 
duced by the neutron branch of the modified Urea process 
and by the nucleon direct Urea process, respectively. 



Recently we (Haensel et al. 200C, hereafter Paper I) 
have considered the bulk viscosity due to the direct Urea 
processes in the npefi matter with superfluid neutrons and 
protons. In the present paper we analyze the bulk viscosity 
of this matter produced by the modified Urea processes. 

Let us remind that damping of stellar pulsations is pro- 
duced also by s hear viscosity (calculated for instance by 
Flowers & Itoh 1976 ). Contrary to the bulk viscosity the 
shear viscosity grows with decreasing temperature T and 
governs the damping of pulsations of cold neutron stars. 
Estimates show (e.g., Paper I) that the bulk viscosity de- 
termines the damping of pulsations in nonsuperfluid stars 
at T ^ 10^ K if direct Urea processes are open in the 
stellar cores or at T ^ 10® K if direct Urea processes are 
forbidden. 

2. Bulk viscosity in non-superfluid matter 

General expression for the bulk viscosity C of nonsuper- 
fluid npe/x matter produced by nonequilibrium (direct or 
modified) Urea processes was obtained in Paper I. We re- 
strict ourselves to neutrino transparent neutron star cores 
(T <^ 10^° K if direct Urea processes are forbidden). We 
will focus on the viscosity which describes damping of stel- 
lar pulsations with frequency ui ~ 10^-10'' s~^. The pul- 
sation frequency is typically much higher than the rates 
of Urea processes. In this high-frequency limit, can be 
written as a sum of the partial bulk viscosities associated 
with all Urea processes. For the modified Urea processes 
{Nl) in npefj, matter. 



c 



Nl 



\Xni\ 



dP 



dXi 



where C,ni is a partial bulk viscosity, P is the pressure, 
— TLi/n\, is an electron or muon fraction (n; and n;, — 
Tin + np being the lepton and baryon number densities, 
respectively), rji — ^J.n—^J'p—^J'h Mi is the chemical potential 
of particle species j, and Xni determines the difference of 
the rates of the direct and inverse reactions of an Urea 
process, Eq. (||): Tni - f^v/ = -XniVi- As in Paper I, 
we assume that in the absence of pulsations the matter 
is in the state of chemical equilibrium [with respect to 
the beta and muon reactions (||)], in which rji — 0. The 
pulsations are supposed to violate this state only slightly, 
\r]\ <c; T. The derivatives in Eq. (^ are taken at chemical 
equilibrium with n^, X^^ and X^^ as independent variables. 
In our notations, dP/dXi and Xni are negative. 

As shown in Paper I, a partial bulk viscosity can be 
rewritten as 



Nl\ 



Cf. Ci 



rib 



drji 

dn.h 



(4) 



In Paper I the factors Cg and were calculated and fitted 
by simple formulae for two model equations of state of 
the npe/i matter proposed by Prakash et al. (1988). These 
equations of state correspond to two different forms of the 
symmetry energy of nucleon matter, S{nh) — {13 [w^/^ _ 
F{u)] + 30F(u)} MeV, where u = rib/no and no = 0.16 
fm^3 jg ^ijg saturation nucleon number density. One has 
F{u) = u for model I; and F{u) = 2u'^ /(I -\- u) for model 
II. 

Thus the problem is reduced to evaluation of the re- 
action rates Tni and the factors Xni- For a modified Urea 
process Nl {h ^ c = ItB = I) 



n 



\ (27r)3j2£i(2^)3 2£,(27r)3 
x(l-/jv')(l-/p)(l-//)/«/iv(27r)4 

x5{Ei-E,)5{Pf-P,)-Y,\M\\ 



(5) 



where Pi is the nucleon momentum (i = n,p, N, N'), pi 
and £; are the electron or muon momentum and energy, 
while Pi, and £y are the neutrino momentum and en- 
ergy. The delta functions S{Ef — Ei) and 5{P f — Pi) 
describe conservation of the energy E and the momen- 
tum P of the particles in initial and final states, i and 
/; |Afp is the squared matrix element of the reaction, 
s = 2 is the symmetry factor which excludes double 
counting of the same collisions of identical particles, and 
fi = {1 + cxp[(£i — fj.i)/T]}^^ is the Fermi-Dirac function. 

To proceed further we need the matrix element M 
which depends on the model of nucleon-nucleon inter- 
action. We employ M calculated by Friman & Maxwell 
( |1979D using the one-pion-exchange model to describe the 
long-range part of the AW interaction and the Landau 
Fermi- liquid theory to describe the short-range part. Their 
result [their Eq. (39)] is that the squared matrix element 
summed over spin states of all particles and averaged over 
orientations of the emitted neutrinos can be written in the 
form: 



J2 |A/|'- 256 Glgl 



r 



ei 



olnPn- 



(6) 



Here, Gp = 1.436 x lO"'" erg cm'^ 
pling constant, gA_ = 1.26 is the axial-vector constant of 
weak hadron current, /'^ w 1 is the ttW- interaction con- 
stant in the p-state in the one-pion-exchange model po- 
tential, and is the pion mass (tt^). Furthermore, un 
describes momentum dependence of the squared matrix 
element in the Born approximation and Pn contains var- 
ious corrections. In principle, these factors are different 
for the neutron and proton branches of the modified Urea 
process (since ri„ ^ rip). Nevertheless, bearing in mind 
the uncertainties of the employed model, we set a„ = ap 
and (3n = /3p. Friman & Maxwell (1979) adopted the val- 
ues a„ = 1.13 (calculated dX p — po) and /3„ = 0.68 (to 
account for correlation effects). In our approximation the 
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squared matrix element in Eq. can be regarded as con- 
stant and can be taken out of the integral. 

Further evaluation of Fat/ is standard (e.g., Shapiro & 



Teukolsky 1983). It is based on phase-space decomposi- 
tion, neglecting neutrino momentum in the momentum- 
conserving delta function and replacing the momenta of 
other particles pi by their Fermi momenta ppi (i = n,p, I) 
whenever possible. Introducing the dimensionless quanti- 
ties 



Mi 



we can present the reaction rate in the form 
F - r(") T- 

>- Nl — >- Nl '^^ 

dxy j dxn dxp dxN dx^i dxi 

X .f{Xn) f{Xp) f{xN) f{xN') f{xi) 
X 5{Xn + Xp + Xn + Xn' + Xi ~ X^ 



-0, 



(7) 



(8) 



(9) 



where rj^| is some typical reaction rate. We have trans- 
formed all the blocking factors (1 — /(x)) into the Fermi- 
Dirac functions f{x) by replacing x — > ~x. The dimen- 
sionless energy integral / for the n and p branches of the 
process is the same. For the reactions involving electrons, 
we have 



7r9 h^^ c8 



3.44 X 10^ 



m^)^ mi I Up >|^^^^7 

m 



p(o) ^ p(o) 

pe ne \ 



■m*\ (3 PFp + PFe - PFn ) 



8pFp PFe 



e 



pe ■ 



(10) 



(11) 



Here, m* and m* are the neutron and proton effective 
masses, respectively, and the step function Ope takes into 
account the threshold character of the proton reaction 
branch: it is allowed {Ope = 1) for ppn < (3pFp + PFe) 
and forbidden {Qpe = 0) otherwise (Yakovlev & Levenfish 
1995). The latter condition holds almost everywhere in 



the core; it may break at p ^ po only for the equations of 
state with small symmetry energy. Notice that Eqs. (^0|) 
and (|ll| ) have to be slightly modified if the direct Urea 
process is allowed (due to modification of the angular in- 
tegrals; e.g., Yakovlev & Levenfish 1995). However these 
modifications are of no practical interest: if allowed, the 
direct Urea process fully dominates over the modified Urea 
processes. 

The rates of the neutron and proton branches of the 
muon modified Urea process are similar. They are ob- 
tained from Eqs. (^) - (|l^) in two steps. First, by introduc- 
ing the extra factor v-p^i/c = (n^/ue)^^^, where v^p. is the 
Fermi velocity of muons. Second, by replacing 8pe — > ©p/i 
that corresponds to the replacement p^e — > PFp- The 
muon step function Op^ allows the proton branch of the 
process at p^n < (3pFp +PFp)- 

The rate F ^v/ of an inverse modified Urea reaction (lep- 
ton capture) differs from the rate of the direct reaction Fjv; 



only by the argument of delta function in the expression 
for /; one should replace ^ — ^ there. The difference of 
the direct and inverse reaction rates can be written as 



Fat; - Fat; = -Ajv; m = 

AI = / dx, xl [J{x, -0~ J{xu + 0] , 



(12) 
(13) 



J{x) = / dxn dxp dxM dxN' dxi /(a;„) f{xp) /(xjv) 



X f{xN') f{xi) 5{Xn+Xp + XM + XN'+Xi ~ x) . (14) 

The function J{x) is evaluated analytically (e.g., Shapiro 
& Teukolsky 1983): 



J{x) 



x^ + IOtt^ x^ + 971" 
24 (e^ + 1) 



(15) 



Thus, the difference ^ ni ~^ Ni in nonsuperfluid matter is 
determined by the only parameter ^ = ^ijT . Furthermore, 
the integral AJ is taken analytically: 



T{i) = 1 



1512 

189g^ 
^ 36 77r2 

Combined with Eq. 

I \ I ^ Nl 
\^Nl\ - —■ 



21C 

3677r4 ^ 18357r6 ' 
|12|), this relation gives: 



(16) 



(17) 



As mentioned above, we do not consider large deviations 
from the chemical equilibrium and restrict ourselves to 
the case \rii\ <C T in which « 1. It is instructive to es- 
timate characteristic times of Urea processes. Taking into 
account that -he — TNe^^Ne ~ A ?7, we can introduce, for 
instance, the time t„ ~ dne/{XeT) required for an elec- 
tron Urea process to shift the chemical potential difference 
by r/ ~ T, within the thermal width of the Fermi levels. In 
this case the electron number density will be changed by 
Sn^ ~ rie (T/fie)- Adopting typical parameters of neutron 
star matter we have ^ 20 T^^ s for direct Urea pro- 
cess (Paper I) and ~ Tg^ months for modified Urea 
processes. Therefore, in any case is much longer than 
typical periods of neutron star pulsations (^ 10^-^ s) which 
justifies validity of the high-frequency approximation used 
in our analysis (see above). 

Finally, combining Eqs. (^) and ( p7| ) we obtain the 
partial bulk viscosity of the npefi matter, Cat/ = Cnio (the 
subscript '0' refers to nonsuperfluid case), produced by 
the modified Urea processes for \rii\/T ^ 1. In standard 
physical units 



CneO 



1512 7r3ai° c8w2 



^ ) (ksTf an f3n 



1.49 X 10 



19 









( ^ 








/ mp 





1/3 



a 



100 MeV 



SpeO 



2 



g cm ^s ^ 



^ , ^\ (3 PFp + PFe -PFnY 
\mt,/ 8 PFp PFe 



e 



pej 



(18) 
(19) 
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CpmO 



,neO 



( Plii 

\PFe 

Cefn* 



(3pFp+PF^-PFr; 



PFeJ 



(20) 

(21) 



where W4 = w/(104 s^^) and Tg = 7/(10^ K). 

For the two equations of state I and II (see above), the 
density dependence of the partial bulk viscosities is fitted 
by simple expressions presented in the Appendix. Figure 
1 shows the total bulk viscosity of nonsuperfluid matter 
versus the baryon number density rih for both models. The 
contributions of the modified and direct Urea processes 
are included; the latter contribution is calculated using 
the results of Paper I. The dotted curves show the bulk 
viscosity for the simplified equations of state based on the 
same nuclear energy but with the appearance of niuons 
artificially forbidden. 

The viscosities for models I and II are similar. The 
muon thresholds are rib = 0.150 fm^'^ and 0.152 fm~^ for 
models I and II. At lower densities the muons are absent 
and the bulk viscosity is produced by the electron mod- 
ified Urea processes (and the dotted curves are identical 
to the solid ones). At densities above the thresholds the 
muon modified Urea processes are switched on and am- 
plify the bulk viscosity. Since wp/^, the muon bulk 
viscosity is switched on without any jump. On the con- 
trary, the bulk viscosity due to the direct Urea processes 
is switched on in a jump-like manner at their threshold 
densities. The thresholds for the electron direct Urea pro- 
cess are rif, — 0.414 fm~'^ and 0.302 fm~^, for models I 
and II. As explained in Paper I, the presence of muons 
lowers the threshold density, mainly due to increasing 
the number density and Fermi momenta of protons. On 
the other hand, the muons lower the bulk viscosity pro- 
duced by the electron direct Urea processes (by decreasing 
He). At larger rib the total bulk viscosity jumps again, at 
rib = 0.503 fm-3 and 0.358 fm'^ for models I and II. This 
time the jump is associated with switching on the muon 
direct Urea process. The contribution of the muon direct 
Urea into the bulk viscosity is even larger than the con- 
tribution of the electron direct Urea (Paper I). With this 
contribution the bulk viscosity becomes larger than in the 
npe matter. 

3. Bulk viscosity of superfluid matter 

3.1. Superfluid reduction factors 

Now we turn to superfluid reduction of the bulk viscos- 
ity produced by the modified Urea processes. As in Paper 
I, we restrict ourselves to the motions of the npen fluid 
in which all components move with same macroscopic ve- 
locity. This assumption reduces equations of stellar pul- 
sations to those of one-fluid hydrodynamics, with a sin- 
gle coefficient of the bulk viscosity; our bulk viscosity ( 
is identical to ^2 of Landau and Lifshitz (1987). The su- 
perfluid reduction of the bulk viscosity is similar to the 
superfluid reduction of the neutrino emissivity. Since the 
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Fig. 1. The bulk viscosity ( for models I and II of nonsuper- 
fluid npefi matter (solid lines), produced by direct and modi- 
fied Urea processes involving electrons and muons, versus the 
baryon number density nt, for T — 10^ K and lo = 10'' s~^. 
Dotted lines show the same bulk viscosity but for the models 
of matter without muons. Jumps of the curves are associated 
with opening direct Urea processes (see the text). 



latter problem has been the subject of extensive studies 



(e.g., Yakovlev et al. 1999) we omit technical details. 

We adopt the traditional assumption (e.g., Yakovlev 
et al. 1999D that the protons form Cooper pairs due to 
smg let-state (^Sq) pairing (superfluidity of type A, in no- 
tations of Paper I). As for the neutrons, they would be 
thought to undergo either the singlet-state pairing or the 
triplet-state ('^P2) pairing with zero projection of the to- 
tal angular momentum of the Cooper pair on the quan- 
tization axis (superfluidity of type B). Numerous micro- 
scopic calculations (see Yakovlev et al. 1999| for references) 
show that the singlet-state pairing of neutrons takes place 
at p ^ po, while their triplet-state pairing is realized at 
higher densities. The critical temperatures Ten and Tcp of 
the neutron and proton superfluids are very model depen- 
dent, and we treat them as free parameters. 

Microscopically, superfluidity introduces a gap S into 
momentum dependence of the nucleon energy, e{p). Near 
the Fermi surface (|p — Pp| Pp), this dependence can be 
written as 



e = p - yJS'^ +vI{p-pyY dX p< pp, 
£ = p + ^J5'^ + Vp{p -pf)^ at p > Pp. 



(22) 



In case A the gap is isotropic (independent of orienta- 
tion of p), (5a — Aa(T), while in case B the gap de- 
pends on the angle •& between p and the quantization axis, 
(5b = AB(r) (1 -I- 3 cos^ I?) 1/2. In both cases A(r) is the 
amplitude which describes temperature dependence of the 
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gap. The amplitude is derived from the standard equation 
of the Bardeen - Cooper - Schrieffcr theory. 

For further analysis we introduce the dimensionless 
quantities (fee = 1 as before): 

A(T) T 5 

11 = _ T — 



y 



(23) 



The dimensionless gap y can be presented in the form: 
UA — Ub — vb -y/l + 3 cos^ The dimensionless gap 
amplitude v depends only on r = T/Tc. This dependence 
is accurately fitted as (Levenfish & Yakovlev |1994 ): 

VA = (1.456- 0.157/ V^+ 1.764/t), 

VB = (0.7893 + 1.188/r). (24) 



Quite generally, the high-frequency bulk viscosity can 
be presented in the form 



Nl 



(25) 



where ^nio is a partial bulk viscosity of nonsuperfluid mat- 
ter [given by Eqs. (|l|)-(|l|)] and = is a factor 
which describes its superfiuid reduction. The factors for 
the modified Urea processes Ne and Nfi are equal since 
the superfluidity affects only nucleons. In non-superfluid 
matter matter, we have R^-^^ = 1 and reproduce the re- 
sults of Sect. 2. 

Thus the problem consists in calculating the reduction 
factors R^^\ Each factor depends on two parameters, w„ 
and Vp, dimensionless gap amplitudes of the neutrons and 
protons. As in Paper I we may assume that the superfluid- 
ity affects noticeably only the factor A^v; in the expression 
for the bulk viscosity, Eq. (^). By making the standard 
angular-energy decomposition in the expressions for the 



reaction rates Tni and r„; 
contains the factor 

5 



one can easily show that A 



Nl 



J{Vp,Vn) = 47r 




(26) 



where dQi is the solid angle element in the direction of 
Pi and A/ is the energy integral given by Eq. (|l3|). This 
integral is the main quantity affected by the superfluidity. 
At ^ <C 1 the integrand of Eq. ( [T3| ) contains the function 
J{x„ — £,) ~ J{xu + ~ dJ(xy)ldxy, where J(xy) is 
given by Eq. (|lj). Thus, for small deviations from equi- 
librium one can transform Eq. (03) to: 



A/ = 4^ / AXyXy 

Jo 



n 



+ 00 



4i: 



(27) 



Here, i — n, p, N, N' , I. 

Let us now label the quantities in the nonsuperfluid 
case by the subscript '0'. The angular-energy decomposi- 
tion yields 

367 TT^^ 



1512 



(28) 



In this case Anio is the angular integral containing the mo- 
mentum conserving delta function, in which the neutrino 
momentum is neglected and the momenta of all other par- 
ticles are placed at their Fermi surfaces. For the processes 
of our study: 



A 



plO 



(47r)^ {3pFp+PFi -pFnf 



8pFpPFl 



(29) 



Notice that Yakovlev & Levenfish ( 1995| ) inaccurately de- 
termined ApiQ. Here we present the corrected expression. 

Generalization of Ajv; to the superfluid case can be 
achieved by introducing the neutron and proton energy 
gaps into Eq. (p6|). Then 



Xm = XNioR'^^\ R'^^Hvp,Vn)-^ 



J 



Anio A/n 



(30) 



where is given by Eq. (26) with 



AI{yi,y2) = dxi, Xu 



n 

.1=1 



dxt f{zi) 




(31) 



Zi = {ei-fi)/T ^ sign(a;i)yif+yf for nucleons, and zi = 
xi . Equation (|3^) is the general expression for calculating 
the reduction factor R'^^^ . 



3.2. Superfluidity of neutrons or protons 

First, let the neutrons be normal while the protons un- 
dergo Cooper pairing of type A. Since the energy gap is 
isotropic, the angular integral in Eq. (|26| ) is separated from 
the energy one, being the same as in the nonsuperfluid 
case, Eq. (p9|). Then in Eq. (|6|) we can put ?/„ — and 
Vp — Vp. Accordingly, the reduction factor depends on the 
only parameter v = Vp. For a strong proton superfluidity 
(r — T/Tcp 1, V 1), the asymptotes of the reduc- 
tion factors for the neutron and proton branches of the 
modified Urea process are: 



r: 



(n) 
pA 



126\/2 



1835 7r5-5 



v^-^ e- 



0.004061 

_5.5 



■ exp 



1.764 



r: 



(p) ^ 6048 
3677r6 



7 



51 



0.1604 



3.528 



16 > 2lA/3-yln(\/3 + 2) 



0.548. 



,(32) 



(33) 



The proton branch is affected by the proton superffuidity 
stronger than the neutron branch. This is natural (e.g., 
Yakovlev et al. 1999) because in the proton branch we have 
three protons which belong to the superfluid component 
of matter, while in the neutron branch there is only one 
such particle. Under the number of reacting nucleons we 
mean their total number (in the initial plus flnal states) 
in the direct or inverse reaction of an Urea process. 
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We have calculated 



and numerically in a 



wide range of v — Vp and proposed the fits which repro- 
duce the numerical results (with the mean errors of < 1%) 
and the asymptotes (p3) and ([3^): 



R 



(n) 
■pA 



,5.5 



exp ( 3.245 - a/(3. 245)2 + t;2 ^ (34) 



a ^ 0.1863 + (0.8137)2 + (0.1310 1;)^ 
b = 0.1863 + V (0.8137)2 + (0.1437 w)2 



r'^^I = exp ( 5.033 - 



V(5.033)2 + {2vy 



(35) 



0.3034 + V (0.6966)2 _^ (0.1437 ^;)2, 



Now consider normal protons and superfluid neutrons 
(superfluidity of type B). Then in Eq. ( p6| ) we may put 
Zp = Xp and the reduction factor depends only on v — Vn- 
For a strong neutron superfluidity (r = T/Tcn ^ 1, w ^ 
1) the asymptotes of the reduction factors are 



R 



in) 
riB 



6048 



R, 



ip) 



367^6 3^ 
42\/3 



0.00720 



1835 



5 ~ 
V e 



3.066x10- 



■ exp 



exp 



2.376 

r 

1.188 



,(36) 
(37) 



In this case the neutron branch is affected by the super- 
fluidity stronger than the proton branch. 

The calculation of R^Bi^n) for intermediate values w„ 
is difficult since there are three superfiuid neutrons with 
anisotropic gaps. An approximate expression for Rl^{vn) 
will be proposed in Sect. 3.3. As for the factor R^^, we 
have calculated it numerically for a wide range of v and 
obtained the fit which reproduces the numerical results 
(with mean error ^ 1%) and the asymptotes: 



d(p) _ 



■exp (2.110- v/(2. 110)2 + u2 



a = 0.1973 + V(0.8027)2 + (0.1257 w)2, 
b = 0.1973 + v/(0.8027)2 + (0.1428 w)2. 



(38) 



The factors i?*^^' which describe reduction of the bulk 
viscosity of npe^ matter by superfluidity of neutrons or 
protons are plotted in Figs. 2 and 3. For comparison, we 
present also the reduction factors R^^^ for the partial bulk 
viscosity produced by the direct Urea processes (Paper I) . 
Figure 2 displays the reduction factors versus T/Tc, while 
Fig. 3 presents them versus v. As seen from these figures, 
the reduction factors may be separated into two groups. 
The first group contains the factors R^^!^ , , R^ , r[^ 



which describe reduction of the reactions with one super- 
fluid nucleon. The other group contains the factors R^jl 

and R^ (asymptote) for the reactions with three super- 
fluid nucleons at once. Notice that the curves inside each 
group in Fig. 3 are much closer to one another than in Fig. 
2. The same situation takes place for the redu ction factors 
of neutrino emissivities (e.g., Yakovlev et al. 1999). 



1 - 



R 



"1 — I — I — i~rT 




direct Urea 
modified Urea 



T/T, 

Fig. 2. Factors R which describe reduction of partial bulk vis- 
cosities of npefj. matter by superfluidity of neutrons or protons 
versus T/Tc- Solid and dotted lines refer to the modified and 
direct Urea processes, respectively. Curves 1 are for the vis- 
cosities due to the neutron modified Urea and the direct Urea 
processes and the proton superfiuidity. Curves 2 are for the 
proton modified Urea and the direct Urea viscosities and the 
neutron superfiuidity. Curve 3 is for the proton modified Urea 
process and the neutron superfiuidity. Curve 4 is the asymptote 
for the neutron modified Urea and the neutron superfiuidity. 

3.3. Superfluidity of neutrons and protons 

If the neutrons and protons are superfiuid at once, cal- 
culation of the reduction factors is complicated. If how- 
ever we do not need very accurate results, we can use 
approximate similarity relations analogous to those sug- 



gested by Yakovlev & Levenfish (1995) and described also 
by Yakovlev et al. ( |1999D for the reduction of neutrino 
emissivities in different reactions: 



^Ab(^P' 



rS!{vp) 



(39) 
(40) 



Here, R'^ , R^^ and i?^^ are the factors which define 
reduction of the direct-Urca bulk viscosity by superfluidity 
of protons, neutrons and protons+neutrons, respectively. 
These factors were obtained in Paper I. Notice that Eq. 

(39) becomes exact for normal protons {vp = 0), while Eq. 

(40) is exact for normal neutrons. We may assume also 
that the reduction of the A^-branch of the modified-Urca 
viscosity by neutron superfluidity can be approximated 
(for V ^ 10) by the factor: R^^2 - R^aM- 

Finally, Fig. 4 illustrates reduction of the bulk viscos- 
ity of npen matter (due to the modified Urea processes) 
with decreasing temperature by superfiuidity of neutrons 
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Fig. 3. Same as in Fig. 2 but versus dimensionless energy gap 
parameter v. 



of type B or protons of type A for rif, = 2nQ and uj = 10^ 
s~^. Thick solid line shows the viscosity of non-superfluid 
matter (cf. with Fig. 1). Thin dashed lines present the 
bulk viscosity suppressed by the proton superfluidity at 
several selected critical temperatures Tcp indicated near 
the curves. The dot- and- dashed line shows the effect of 
neutron superfluidity (Ten = 10^*^ K) for normal protons. 
We see that superfluid reduction of the bulk viscosity de- 
pends on temperature, superfluidity type, and the critical 
temperatures Ten and Tcp. If Ten and Tcp arc not higher 
than 10^" K for rn, ~ 2 no, then the superfluid reduc- 
tion cannot be very large, say, for T ^ 10^ K. It cannot 
reach more than two orders of magnitude in the case of 
superfluidity of one nucleon species or six orders of mag- 
nitude if n and p are superfluid at once for T = 10^ K 
at Ten = Tcp = 10^° K (Fig. 5 of Paper I). The reduction 
grows exponentially with further decrease of T. 



4. Conclusions 

We have calculated the bulk viscosity due to non- 
equilibrium modified Urea processes in neutron star cores 
composed of neutrons, protons, electrons, and muons 
{npeji). We have considered non-superfluid matter (Sect. 
2) and described also reduction of the bulk viscosity by 
superfluidity of neutrons and protons (Sect. 3). In combi- 
nation with our previous work (Paper I) we have obtained 
realistic description of the bulk viscosity provided by the 
modified and direct Urea processes in the npe/i matter 
of the neutron star cores. The results can be used for 
studying damping of pulsations in neutron stars and grav- 
itational radiation driven instabilities in rotating neutron 
stars. Strong superfluidity of neutrons and protons reduces 
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Fig. 4. Bulk viscosity of superfluid npeji matter (model I) 

produced by the electron and muon modified Urea processes 
at the baryon number density rib = 2 no and uj — 10'' s^' as 
a function of T. Thick solid lino corresponds to non-superfluid 
matter, dashed lines axe for matter with superfluid protons 
(Tcp = 10^°, 10^-® and 10^'^ K) and normal neutrons, while dot- 
and-dashed lino is for matter with superfluid neutrons {Ten = 
10^^° K) and normal protons. Curves are labeled by Ig Tcp or 

IgTcn. 

the bulk viscosity and creates favorable conditions for the 
development of these instabilities. 
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Appendix 

To calculate the bulk viscosity from the equations ob- 
tained above one needs to know the factors Cg and 
and the number densities n„, rip, rxe, n^, and rib. Other 
functions of density such as Tj^i are expressed through 
the Fermi-momenta of various particles ppi = ?i(37r-^nj)^/^ 
and hence through rij. The practical expression for cal- 
culating Ce and Cp for any given equation of state was 
obtained in Paper I [Eq. (19)]. The number densities rij 
are usually available for a given equation of state. It is 
sufficient to specify rip, and Ub, since nn = rib — np and 
rig — rip ri^. 

For example, consider two model equations of state I 
and II used in Paper I and in the present paper. The fits of 
the quantities qi = (ni/no)"^/^(C;/100MeV)^, I = e and /x, 
as functions of the baryon number density rib are given by 
Eq. (36) of Paper I. They yield practical expressions for 
Ci{nb) provided the particle number densities are known. 
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Thus, we must fit the number densities Ui ~ ni{nb). We 
present the proton number density as Up — UpQ + Aup, 
where UpQ is the proton number density for the case in 
which creation of muons is artificially forbidden. For the 
equations of state I and II the muons appear at rib = n-bp, = 
0.15000 and 0.5122831 fm'^, respectively. Therefore, we 
have Up = Upo = Ue and = for Ub < rib^. 

We have fitted UpQ in the density range from Ub — 0.048 
fm^^ = 0.3 no to 1.92 fm^^ ^ 22no as 

npo = r] + an^^/{l+-fn^), (Al) 

and we have also fitted Aup and in the range from 
Hb — nbfj, to 1.92 fm~^ as 

Aup = aSn^-'{l + dSn'')/{l + b5n''), 

= ASn^-^{l + DSn'^)/{l + BSn^). (A2) 

Here, n = Ub, Sn ^ Ub — rib^,, and all number densities are 
expressed in fm^^. The fit parameters are a = 0.7295, p = 
2.503, 7 = 2.604, S = 1.307, r/ = 0, a = 0.1306, b = 1.887, 
c = 1.509, A = 0.1614, B = 0.1785, C = 1.235, d = D = 
for model I; a = 3.64, f3 = 3.374, 7 ^ 12.07, d = 2.096, 
77 = 2.673 X 10-5, a = 0.1919, b = 4.946, c = 1.279, d = 
0.1558, e = 0.9253, A = 0.2376, B = 3.294, C = 1.468, 
D — 2.194, E = 1.201 for model II. The maximum errors 
of the fits of Hp and do not exceed 1.3% in the indicated 
density ranges. The presented fits give all functions versus 
rifc. The fit expressions of nb via mass density p for three 
versions of each equation of state I and II corresponding 
to three different values of the compression modulus of 
saturated nuclear matter, K = 120, 180 and 240 MeV, 
are given by Eq. (37) of Paper I. 
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